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ABSTRACT 

We present a comprehensive study of the evolution of Luminous Red Galaxies (LRGs) 
in the latest and final spectroscopic data release of the Sloan Digital Sky Survey. We 
test the scenario of passive evolution of LRGs in 0.15 < z < 0.5, by looking at the 
evolution of the number and luminosity density of LRGs, as well as of their clustering. 
A new weighting scheme is introduced that allows us to keep a large number of galaxies 
in our sample and put stringent constraints on the growth and merging allowed by the 
data as a function of galaxy luminosity. Introducing additional luminosity-dependent 
weighting for our clustering analysis allows us to additionally constrain the nature 
of the mergers. We find that, in the redshift range probed, the population of LRGs 
grows in luminosity by 1.5-6 % Gyr -1 depending on their luminosity. This growth is 
predominantly happening in objects that reside in the lowest-mass haloes probed by 
this study, and cannot be explained by satellite accretion into massive LRGs, nor by 
LRG-LRG merging. We find that the evolution of the brightest objects (with a K+e- 
corrected M r .o.i ^ —22.8) is consistent with that expected from passive evolution. 
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1 INTRODUCTION 

Luminous Red Galaxies (LRGs) are useful probes of large 
scale structure. They are bright and have the potential 
to be easily separated by colour, allowing a uniform sam- 
ple to be selected with which to map large volumes. They 
are also strongly biased with respect to the matter dis- 
tribution, making them particularly attractive for stud- 
ies that aim to use Baryon Acoustic Oscillations (BAO) 
in the ga laxy distribution to study the expansion of the 
Universe dSeo fc Eisensteinll2003l:|Blake fe GlazebrookllioO^: 



iHu fe Haimanl 120031 ; iMatsubaral I2004T ). Because of this, 
LRGs are being targeted by the Baryon Osc illation Spec- 
troscopic Survey (BOSS; ISchlegel et "all I2009T ). part of the 
SDSS-III project, in order to map the expansion rate out to 
z < 0.7. 

In addition to allowing measurements of cosmologi- 
cal expansion by using the BAO signal as a standard 
ruler, galaxy surveys also measure the gro wth of struc- 
ture through redshift -spa ce distortions ( e.g., iKaiserl [l987l ; 
iDavis fc Peebles|[i"983l and lHamiltonll 1998 for a review). Al- 
though current BAO and z-space distortion analyses are 
generally robust to an evolving galaxy bias, f uture experi- 
ments will be more sensitive to these effects fe.g. lSmith et al] 
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l2007h . Observing a single population of passively evolving 
galaxies offers an attractive way of minimising bias effects. 
In particular, for such a sample, the redshift evolution of 
bias and its effects on cosmological probes can be more easily 
modelled. Uniformity in terms of selection is also advanta- 
geous for cosmological studies of large-scale structure, which 
are often simplified if the relationship between the galaxies 
and the underlying density field is uniform across the sam- 
ple. It is therefore extremely interesting to ask whether a 
uniform, passively evolving population of galaxies can be 
found. 

LRGs are the most likely candidate for such a sample, 
as they are traditionally assumed to form a single population 
of galaxies, which assembled at high-redshift and has been 
passively evolving since. As red galaxies dominate the stel- 
lar mass budget in the Universe, understanding how they 
assemble and grow is a key question in theories of galaxy 
formation and evolution. Driven both by galaxy evolution 
and observational cosmology, passive evolution of LRGs has 
been tested extensively. Traditionally, this is done by looking 
at three main observables and their evolution with redshift: 
the total comoving number density, the luminosity function 
and density, and the large or small-scale clustering. 

Firstly and most simply, passive evolution predicts that 
the number density of objects must be conserved as a func- 
tion of time. This has been tested bv lWake et all (|2006l ) , who 
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found that the number density of LRGs brighter than a given 
luminosity is conserved between z = 0.55 and z = 0.2, if one 
accounts for the passive fading of the stellar populations. 

Secondly, we can go further and look at the evolution of 
the lu minosity density as a function of redshift. Wake et al.l 
(2006) used the 2SLAQ sample to investigate whether the 
luminosity function (LF) at z = 0.2 and z = 0.55 is consis- 
tent with a purely passive evolution. They found that this 
is indeed the case. Furthermore, they also looked at the ef- 
fect of major LRG-LRG mergers in the LF. They found that 
the observed LFs are ful ly consistent with a merger-free his- 
tory. [Brown^etj]I] l|2007l ). focusing on the B-band luminosity 
density, have found a slight departure from pure passive evo- 
lution, and that around 80% of the stellar mass in massive 
red galaxies (> AL*) today was already in place by z = 0.7. 
This results in a modest growth of w 3% Gyr -1 in mass 



or luminosity. ICool et "all (|2008l ) put a similar limit on the 
growth of very bright red galaxies (> 3L*), and showed that 
bright red galaxies can not have grown by more than 50% 
since z = 0.9, or w 6.8% Gyr "\ 

These studies face a fundamental problem that, while 
we can match numbers of objects and the luminosity func- 
tion to those expected for a particular evolutionary theory, 
we cannot test whether this is a coincidence. The match is 
necessary to test the theory, but is not sufficient to show that 
the theory must be true. In particular, for passive evolution, 
one can potentially have galaxies entering and/or leaving the 
red sequence at any point during the tested timeline, and 
several effects can change the number or luminosity density 
of a sample. This can happen by: 

(i) blue galaxies quenching star formation and becoming 
red; 

(ii) fainter red galaxies merging together without trigger- 
ing star formation (SF) to enter the sample at lower redshift; 

(iii) LRGs merging with other LRGs to decrease the num- 
ber density at low redshift 

(iv) LRGs leaving the red-sequence due to merger- 
induced SF. 

The lack of massive blue galaxies makes i) and iv) weak 
propositions. Furthermore, LRG-LRG mergers tend to be 
dry, given that massive red galaxies have little in the way of 
gas and SF. But ii) and iii) almost certainly play a role in 
the overall LRG population evolution. The question is: how 
much so? 

One can hope to distinguish between these scenarios 
by investigating the clustering of the LRGs. The two-point 
correlation function at very small scales can be used to in- 
fer merger rates, assuming all LRG-LRG pairs closer than a 
given distance will eventually merge within a t imescale given 
by th e orbital or dynamical friction time. iMasiedi et al.l 
(2006) put an upper- limit on the LRG-LRG merger rate 
of « 0.6 x 10 4 Gpc 4 Gyr ~ 1 . This very low rate suggests 
that mergers between LRGs i s not the prima r y cau se of evo- 
lution in the population. In IMasiedi et all ( 20081 ) the au- 
thors extended this reasoning to other types of galaxies, by 
computing the cross-correlation between LRGs and both red 
and blue galaxies, of different luminosities. They found that 
most of the luminosity brought to LRGs comes from red 
galaxies via dry mergers and they put an upper limit on 
the growth of LRGs at 1.7± 0.1 h % Gyr -1 at z a 0.25. 
These upper limits sit of the low end of other measurements, 



but it should be pointed out the IMasiedi et all |2006l . l2008t ) 
measurements only take into account LRG growth that in- 
volves an already existent LRG - i.e., it fails to account for 
new LRGs being formed from the mutual mergin g of fainter 
(red) galaxies. Similarly, |Pe Propris et all (|2010l ) looked at 
dynamically close pairs in the 2SLAQ sample (concentrating 
on pairs of galaxies in 0.45 < z < 0.65), and f ound a merger 
rate o f 0.6 x 10 4 Gpc 4 Gyr _1 , consistent with lMasiedi et aU 
(2006) and reaching a similar conclusion. 

Finally, one can look at the evolution of the cluster- 
ing of the LRGs and ask whether it is consisten t with 
pas sive-evolution. At i ts simplest, both I White et all (|2007T l 
and lWake et al" I <|200Sl ) found no evolution in the clustering 
amplitude as a function of redshift. This alone is a signifi- 
cant problem for passive evolution, as a scale-independent 
and deterministic bias predicts evolution in the clustering 
strength. To go further, such analy ses requires a frame - 
work for modelin g this evolution. IConrov et al.l d2007t) ; 
IWhite et al l (|2007l ): lBrown et al l |2008l ): IWake et al.l J200?li 
have all, albeit in slightly different fashions, used the halo 
model to support their interpretation. Consistently, these 
studies find that, within the halo model framework, where 
there are too many satellites at low redshift if one is to 
assume passive evolution from the observed clustering at 
high-redshift (33 to 50 per-cent of the satellites must dis- 
appear). These galaxies must either merge with the central 
galaxy or b e disrupted and beco me part of the intra-cluster 
light (ICL). IConrov et all (|2007l ) suggest that a significant 
amount of the stellar mass in the merged haloes never makes 
it to the central galaxy, and that this is a way to reconcile the 
observed lack of evolution in the mass or luminosity function 
of LRGs with clustering studies that indicate that a large 
number of mergers happened since z=l. This argument is 
given some strength by the fact that they can match the 
total stellar mass budget of clusters, and by the prediction 
of a significant amount of ICL around satellite g alaxies that 
was seen in the Virgo cluster |Mihos et al. I l2005l ). 

Comparison between all of these studies is made very 
difficult by the different selection criteria, redshift range, and 
perhaps most importantly the number density of the sam- 
ple, which in turn determines the luminosity/masses probed. 
Evolution of galaxies in the red sequence is heavily depen- 
dent on luminosity, with lower-luminosi ty galaxies sufferin g 
from significant evolution since z=l (e.g. lBrown et ai1l2007l ). 

1.1 This work 

In this work we present an analysis of the evolution of LRGs 
in the complete sample of galaxies o bserved using the orig i- 
nal SDSS LRG targeting algorithm l|Eisenstein et al.ll2001f ). 
We analyse the number and luminosity density of these 
galaxies, the evolution of the luminosity function, and clus- 
tering strength. 

Our methodology significantly differs from previous 
works in two ways. Firstly, we introduce a new weighting 
scheme that allows us to match galaxies at low- and high- 
redshift whilst keeping most of the galaxies in the sample. 
Secondly, as part of our clustering analysis we analyse a 
luminosity-weighted power-spectrum and its evolution. On 
large-scales, this statistic is insensitive to the merging of 
galaxies within the sample assuming conservation of light, 
and provides information on the causes of evolution and the 
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population of galaxies affected. By comparing clustering for 
luminosity and number-density weighted samples, we bypass 
the need to introduce the halo model in order to analyse how 

evolution is occurring. 

Throughout the paper, we use the [Maraston et al.l 

(2009) fiducial model for the expected colour evolution of 
an LRG, assuming passive evolution (M09). The stellar evo- 
lution is the only model dependence of this work, apart from 
dynamical passive evolution that we aim to test. Our choice 
of model is further discussed in Sections 12.11 and l6l 

Our paper is organised as follows: in Section [2] we 
describe our samples, including selection, K+e-corrections, 
weighting-schemes and sample matching; in Section[3]we ex- 
plain our method for the computation of the galaxy power- 
spectrum and in Section|4]its predicted evolution in the pas- 
sive scenario; results are presented in Section [5] we explic- 
itly consider uncertainties in the stellar model in Section [6] 
we consider the implications for the LRG population and 
present our interpretation in Section [7] a comparison with 
previous work is done in Section [8] and finally we summarize 
in Section [9] 



2 THE LRG SAMPLE 

We analyse the SDSS spectroscopic LRG sample, whose 
selection in colour and luminosity was described in 
lEisenstein et al.l l|200l[ ). The upper panel of Fig.[T]shows the 
distribution of all of the galaxies in terms of colour. The red 
line shows the track predicted by the M09 model for a galaxy 
made up solely by stars that are 12 Gyr old at redshift of 
zero, and that have been left to passively evolve since. The 
lower panel of Fig. [T] shows the redshift distribution of the 
galaxies. In this paper we adopt the following nomenclature: 
SDSS model magnitudes are written as u, g, r, i, z; SDSS 
petrosian magnitudes as up, gp,rp,ip, zp and magnitudes 
given by the fiducial model of M09 as u m , g m ,rm,r m ,i m , z m . 

The M09 model is used to match samples at low and 
high redshift, and in calculating rest-frame absolute magni- 
tudes. The SD SS LRG selection, which was base d on the PE- 
GASE models l)Fioc fc Rocca-Volmerangell 19991 ). is only im- 
portant in that this defines the galaxies whose evolution we 
are testing. I.e. we can decouple the selection of the galaxies 
whose evolution we are testing with the model used to test 
that evolution. Essentially, in t his paper, we test wheth er the 
SDSS LRGs selected using the lEisenstein et all l|200ll) cuts 
are dynamically consistent with passive evolution, assuming 
a model for the colour evoluti on that we know is a ble to 
fit at least part of this sample l|Maraston et al.| [2009') . One 
may also want to explore whether samples of galaxies that 
are further or closer to the fiducial model, in colour-colour 
space, follow passive evolution more or less closely. We leave 
that for a follow-up paper. 



2.1 K+e corrections 

We have used the observed colours and the fiducial model 
to compute the K- and K+e-corrections used throughout 
this paper. The M09 model is the only published model that 
matches the colour evolution of LRGs over our redshift range 
of interest, and therefore it is the only viable option for this 
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Figure 1. Top: Colour-colour plane for LRGs. The black points 
are a sub-set of 10,000 galaxies, and the blue contours are rep- 
resentative of the number density of points (as calculated from 
the full sample). The red line corresponds to the fiducial model of 
Maraston et al (2009). Bottom: Redshift distribution of the orig- 
inal sample. Below a redshift of 0.15 the LRG selection is heavily 
contaminated, and these galaxies are later discarded. 



work. The solution offered in M09 stems from two major ad- 
ditions with respect to the previous literature: the inclusion 
of a sub-dominant metal-poor stellar component (as opposed 
to a yo ung stellar component) , and a change in the stellar 
library. IConrov fc GunrJ (|2009l ) have raised some concerns 
about the latter because they, with t heir SPS method, coul d 
not replicate the changes reported in lMaraston et all p009). 
This inconsistency will be considere d further in Maraston e t 
al. (2010, in prep). Additionally, in iMaraston et all (2009), 
the authors explicitly assume dynamical passive evolution 
of the data to which the model is fitted. If this assumption 
is relaxed one may well find that the inclusion of a metal- 
poor stellar component is no longer the only - or indeed the 
best - solution to the problem. The assumption of passive 
dynamical evolution, however, is the very hypothesis this 
paper aims to test, and therefore perfectly consistent with 
our analysis. Thus use of M09 models, which are the best fit 
to the ridge line of a sub-sample of SDSS selected LRGs, is 
consistent with the test we are performing. 
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The fiducial model provides L\(t age ), the luminosity 
per unit wavelength of an LRG of age t age . To make our 
analysis more straightforward, we K+e-correct all galaxies 
to a common redshift of z c — 0.1, and calculate corrected 
absolute magnitudes as 
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5 logj 



DL{zi) 
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Ke(z, z c ) = 
= -2.51og 10 
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X is in the observed frame and A e in the emitted frame. T\ is 
the SDSS's r-band filter response, and L\{z) the luminosity 
of a galaxy at redshift z, given the fiducial model. Tradition- 
ally, calculating K-corrected rest-frame magnitudes is done 
to the mean redshift of the sample, and is accompanied by 
a change in the filter's wavelength such that, for galaxies at 
that redshift (which should be the majority of the sample) , 
the K-correction is independent of the galaxy's spectrum 
(often not known). However, note that Equation [2] gives a 
fixed K+e correction for a given redshift, and is independent 
of the actual observed colours or spectrum of each galaxy. 
The assumption is that the model is a good description of all 
galaxies in the sample. This makes the choice of the common 
redshift at which to normalise spectra, z c , purely arbitrary. 

A comparison of the rest-frame r-band K-corrected ab- 
solute magnitudes obtained using this model and the ones 
obtained using the code K-correct (jBlanton fc Roweisll2007t ) 
can be seen in Fig. [5] The scatter can be explained by the 
fact that K-correct does not use a fixed template, but rather 
fits a spectrum to the photometry in order to find the K- 
corrections. There is a small offset of around 0.03 magni- 
tudes, which is roughly constant with magnitude, and the 
agreement is generally good. 
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Figure 2. Top: The difference in k-corrected only rest-frame ab- 
solute magnitudes in the r-band when obtained with either k- 
correct (K07) or the M09 model. There is a small offset of around 
0.03 magnitudes, which is roughly constant with magnitude. Bot- 
tom: Distribution of K+e corrected rest-frame magnitudes in the 
r-band for all galaxies at z > 0.15. 



2.2 Sample matching 

In order to analyse the evolution of the LRG population, 
we construct a number of pairs of samples, with one sample 
in each pair at low and one at high redshift. The primary 
difficulty in doing this is selecting a sample at low-redshift 
that matches, in terms of individual galaxy properties, the 
evolved product of the sample at high-redshift. This sample- 
matching process has to take into account four redshift de- 
pendent effects: 

(i) the intrinsic evolution of the colour and brightness of 
an LRG, 

(ii) the varying errors on galaxy colour measurements, 

(iii) the varying photometric errors in any one band. 

(iv) the varying survey selection, 

Our correction for (i) is unashamedly model-dependent. 
We work under the assumption that M09 is a good descrip- 
tion of the colour evolution of the stellar populations, and 
assume passive dynamical evolution for the galaxies in our 
sample - the latter does not pose a problem for this work 
since it is the suitability of this model we are trying to inves- 
tigate. We include an evolving colour scatter term to allow 
for (ii). This is described in Section \2. 2. II 



We show that (iii) is a negligible effect in Section fe.Sl bv 
creating mock samples. 

Traditionally (iv) has been satisfied by removing galax- 
ies that c ould not have been obse rved in the other sample in 
the pair (|Wake et al.ll2006l . [20081 ). However this is wasteful 
in that, in order to fully match samples including distribu- 
tions within each, we need to remove galaxies that could 
not have been observed across all of both samples. In con- 
trast, our approach for (iv) is to construct a set of weights 
that assures that each population of galaxies - in terms of 
colour and absolute magnitude - is given the same weight 
in the high and low redshift samples. For a catalogue whose 
selection is only based on a magnitude limit, the traditional 
Vinax estimator could be used for this instead. We explain 
our weighting scheme in Section [2.2.31 which differs in that 
it is designed to be optimal in the limit of Poisson errors. 

2.2.1 Predicting LRG colours and magnitudes with 
redshift 

In order to match samples, we need to predict the colour 
and magnitude that an observed LRG would have if it was 
moved to a range of redshifts following passive evolution. 



LRG evolution in DR7 5 



0.25 




0.0 0.2 0.4 0.6 

redshift 

Figure 3. o g - r (solid line) and a r —i (dashed line), calculated 
with respect to the fiducial model of Maraston et al (2009), as 
a function of redshift. The red lines are the standard deviations 
obtained after correcting for the fact that we are only sampling 
a finite range in colour and are the ones we use in this work. The 
vertical lines represent the redshift range used. 



This will depend on i) the observed colours and absolute 
magnitude; ii) the fiducial evolutionary model. In order to 
match observed samples we also need to account for the evo- 
lution of the colour-colour scatter that has a clear redshift 
dependence. 

We are assuming a single evolutionary model, so for any 
given redshift the predicted colours are fixed. To calculate 
the colours of an LRG at any redshift, we must include some 
sort of departure from the fiducial model that takes into 
account the observed scatter. 

We start by calculating the scatter with respect to the 
model, as a function of redshift, for each of the two colours. 
This can be seen in Fig. O where we show a g - r (z) and 
a r -i{z). Due to target-selection cuts in the colour-colour 
plane, we are not guaranteed to fully sample the distribution 
of g — r and r — i, which may affect our estimate of a g - r (z) 
and a r -i(z). We are in the regime where this effect is small 
and we assume that the estimation of the mean is unbiased. 
We therefore apply a correction to a given by 



Ao z 



(x — (x)) p(x)dx + 



(x — (x)) p(x)dx (3) 



where a and b represent the limits sampled by the data, 
and p(x) is the probability distribution which we assume 
to be Gaussian. In practice we estimate a and b from the 
data at each redshift (finding that \a,b\ S> u), and use the 
uncorrected value of a to estimate p(x) and evaluate the 
integral. Fig. shows that the correction is small within the 
redshift range probed. 

Having calculated the scatter, we require that the colour 
at a redshift z' of a galaxy observed at zo departs from the 
fiducial model in a way such that the ratio of the distance 
from the data point to the model, A(g — r), and the scatter 
at that 2q is maintained at any other redshift. Explicitly: 



(<? - r)(*0 = {g- r) m {z') + A(g - r)(z )- 



(4) 



&g-r{za) ' 

and similarly for r — i. Note that we are assuming that the 
observed scatter in colour is mostly intrinsic or due to pho- 
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Figure 4. K+e corrected absolute magnitude as a function of 
redshift. 



tometric errors. I.e., that both cuts are wide enough that 
they sample the LRG population almost completely. This is 
not strictly true, but we find that it is indeed almost exactly 
true, as can be seen by the corrections applied to (J g ~ r and 
(Jr-i in Fig. [3] 

We calculate petrosian r-band magnitudes and surface 
brightness as a function of redshift as 

D L (z') 



r-p(z') = Mrfl.i + 5log 1 



lOpc 



fj-so(z') = rj»(«') + 2.5 log 10 [2^50 («')]. 



M*o)§^ with Da being 



+ Ke(z',0.1) (5) 
(6) 

The K+e-corrections are given by Equation [2] Equation [6] 
assumes that the physical size of a galaxy does not change 
with redshift, and $so(z ) 
the angular diameter distance. 

Note that as it is written, the predicted evolution of 
r p (z) does not take into account the changing photomet- 
ric error with redshift. The observed scatter in r p (z) is due 
both to the increase in the photometric errors and the scat- 
ter in intrinsic luminosity. It is hard to disentangle the two 
on a galaxy by galaxy basis, so instead we test the effect 
of assuming a constant error with redshift on our sample 
matching and selection methodologies. The effect is found 
to be negligible, and we give more details in Section \2. 31 

This gives us the colours, petrosian magnitude and sur- 
face brightness as a function of redshift, which allows us 
to know where each galaxy would be observed within the 
survey. Recall that for each galaxy we also have a redshift, 
and a K+e corrected r 01 -band rest-frame absolute magni- 
tude. The redshift - absolute magnitude plane can be seen in 
Fig. 13] This figure clearly shows that the SDSS colour cuts 
result in a sample that doesn't easily lend itself to a volume- 
limited analysis. The two sharp diagonal edges correspond 
to the two rp-band limits in cut I and cut II at rp < 19.2 
and rp < 19.5 respectively. 



2.2.2 Photometric errors 



As mentioned in Section 12.2.11 there is an added subtlety 
in this selection that comes from the fact that photometric 
errors increase with redshift. We do not model this in r p 
when we construct and match our samples. That puts us in a 
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potentially vulnerable position with respect to a Malmquist 
type of bias, in which a slope in the number density with 
luminosity may result in more faint objects being scattered 
into a sample than bright objects are scattered out. Note 
that if the error was constant with redshift our weighting 
scheme would perfectly correct for that because we match 
the weight of each population at high and low redshift, for 
any given magnitude cut - so we would match any excess or 
deficit of objects. 

To check the effect of assuming a constant error we do 
the following. We estimate the true luminosity function of 
LRGs by using a Vmax weight on our full catalogue, and 
use it to construct high and low-redshift mock catalogues in 
the magnitude ranges that we adopt for our final catalogues 
(see Section r2.3|) . We do this twice for each catalogue: once 
accounting for the changing error with redshift in r p , and 
once assuming a constant error. The effect in the number 
and luminosity density in each catalogue between the two 
methods is between 0.05% and 0.3%, depending on the mag- 
nitude cuts. The effect is small enough that we ignore it for 
the rest of this paper. 



2.2.3 Weighting scheme 

In this paper we only consider two redshifts slices, which oc- 
cupy a volume Va (at high-redshift) or Vb (at low-redshift). 
Using the analysis described above, we can calculate the 
maximum volume within each redshift slice that a galaxy 
could have been observed in, given its measured colours and 
redshift. A traditional Vmax correction to match the sam- 
ples at high and low redshift would use this information to 
up-weight galaxies by the reciprocal of the fraction of the 
volume of the respective redshift slice within which a galaxy 
could be observed, i.e. we should apply a weight to each 
galaxy 



Valic 



(7) 



where V ma x,i is the volume the galaxy could have been ob- 
served in a redshift slice, given the evolution of its optical 
properties and the characteristics of the colour cuts, and 
y slice is the total volume of the redshift slice the galaxy falls 
in. Galaxies that can be observed across all of the slice in 
which they fall are given a weight of unity, but galaxies that 
could only have been observed in part of the slice are given a 
larger weight, to compensate for identical galaxies that exist 
at other redshifts but fail to meet the survey selection cri- 
teria. Such a weighting scheme fails to remove any galaxies 
that are not seen at all in the other redshift slice, leaving 
the match between the two slices unbalanced in terms of 
galaxy properties. Furthermore, we also risk massively up- 
weighting populations that only exist in very small numbers 
in a slice, and Poisson errors of such galaxies may dominate. 

Instead, we construct a new scheme that keeps the total 
weight of each galaxy population the same in the different 
redshift slices by down-weighting galaxies based the min- 
imum fractional volume that the galaxy could have been 
observed in for both slices. Explicitly, for a galaxy in Va we 
calculate 
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Figure 5. A schematic representation of our weighting scheme, 
given by Equations [8] and [9] We give three different examples. In 
each case, the grey stripe represents the volume a given type of 
galaxy could be observed in (see main text for details). 



and similarly for a galaxy in Vb- 



V b 



V 1 



V' 



Vi 



V A ' V B 



(9) 



Where the traditional V max estimator would up-weight 
galaxies, we instead down-weight the corresponding galaxies 
with the same properties in the other slice. 

In Fig. [S] we schematically show three representative 
populations of galaxies. Here population refers to a set of 
galaxies that have similar colours and luminosities. A galaxy 
in population 1 that is observed in Va, for example, is given a 
weight of unity, according to equation [8] A traditional Vmax 
weight would up-weight this galaxy by Va/V^^^, but our 
approach instead down-weights the population 1 galaxies 
seen in Vb- This automatically takes care of problematic 
situations like population 2, which is constrained to one of 
the redshift slices - in this case it would be given a weight 
of zero. Population 3 simply represents the galaxies that we 
can see everywhere in the survey - in this case the weight 
is always unity meaning we are of course volume-limited 
for this particular population. The weighting is done on a 
galaxy- by- galaxy basis. 

We have to be careful in the interpretation of results 
where we use our weighting, compared with the traditional 
Vmax weighting. Whereas the latter gives us the means to 
correct for incompleteness and yields true space densities, 
the former should be thought as a weighting scheme rather 
than a completeness correction. This means that weighted 
number and luminosity densities calculated in this way are 
still potentially volume incomplete, but the populations are 
weighted in such a way that they are equally represented 
at both redshifts. We can compare the distribution of to- 
tal weighted luminosity for the two slices, but we cannot 
interpret these functions as giving the luminosity density. 

The advantage of this weighting scheme over simply 
matching samples by removing galaxies that do not match 
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the selection criteria in the other slice is that we keep a large 
proportion of the sample - in fact, the number of galaxies 
thrown out due to the situation of population 2 in Fig. [5] 
is very small. These galaxies contain no information about 
galaxy evolution between the two samples as galaxies of 
that type (luminosity, colour) only exist in one of the sam- 
ples. Consequently, no information is lost by removing these 
galaxies. It is easy to see that the remaining galaxies are 
weighted optimally in terms of Poisson statistics: the impor- 
tance of each sub-population is dependent on the smallest 
number of galaxies of that type in either sample, and is given 
a weight accordingly. We present the final size of our sam- 
ples, after applying the sample selection method described 
in the next Section, in Table [2731 



2.3 Sample selection 

In the rest of this paper we consider only the galaxies with 
redshifts in [0.15, 0.5]. We start by selecting galaxies that are 
brighter than a given chosen absolute K+e corrected magni- 
tude, M'*cu , which can be varied. The sample is split at the 
median redshift, to ensure that we have the same number of 
objects in each slice to begin with. As in Section 12.2.31 we 
have labelled the high-redshift sample as sample A, occupy- 
ing a volume Va and the low-redshift sample as sample B, 
over a volume Vb- 

We then calculate the integrated weighted luminosity of 
the sample as 



Ia 



N(L)LdL 



E 

Li>L mir 



LiWi, 



(10) 



where N(L)dL is the number of LRGs with luminosity in 
[L, L + dL]. In practice we simply use L = 10 _M *" 7/2,5 as 
a proxy for luminosity, which is essentially a luminosity in 



V2.5 



Mlim.A • i • 

r • is arbi- 



arbitrary units and L m in,A = 10~ 
trarily chosen and can be varied so that we may learn about 
the evolution of LRGs as a function of magnitude. We can 
the compute a weighted luminosity density as 



/ Ia 
Va' 



(11) 



This is the luminosity density expected in the r 01 -band, 
at z = 0.1, from LRGs in the high-redshift slice, weighted 
to match the galaxies observed at low redshift. Alterna- 
tively and in an identical fashion, we can also calculate the 
weighted number density of the sample as 



N(L)dL = J2 



and define a total weighted number density as 



n A 



Na 
Va ' 



(12) 



(13) 



We now wish to construct a sample at low redshift that 
matches either the luminosity or number density of the sam- 
ple at high redshift. To do this we find L m in,B such that 



1 f°° 1 
— / N(L)LdL = — V L iWi = £ A 

Vb J Lmin,B Vb "H „ 



or, alternatively, 

1 r°° 1 

^ J N(L)dL = — Wi = nA - ( 15 ) 



TIB 



We can always associate a M^™{ B with L 



min,B , 



but 



note that it will be different from M l *™[ A , even for samples 
matched by luminosity-density. This on itself is a potential 
estimator of how much a sample deviates from passive evo- 
lution. However, in this paper we will instead concentrate on 
weighted number and luminosity densities alone - a compari- 
son between the luminosity-density obtained when matching 
on number-density and vice- versa may immediately tell us if 
the data does not support with passive evolution, we look at 
this in Section r5.ll For the rest of this paper, when we con- 
sider the limiting magnitude of the sample we are referring 
to the magnitude cut at high-redshift. For reference, we give 

M lira, A -i * rlim.F 

r,0.i and M r,0.i 



for each of our samples in Table 12.31 



2.3.1 Choice of za and zb 

The entirety of this work is based on comparing two red- 
shift slices. However, both the boundary between these two 
slices and the redshift distributions within each slice de- 
pend heavily on the faint magnitude chosen. As we include 
fainter galaxies, we predominantly enrich the sample with 
low-redshift galaxies, producing a gradient of galaxy den- 
sity across each sample that will depend on the magnitude 
limit of that sample. The picture is further complicated as 
we weight galaxies as described in Section 12.2.31 We must 
therefore consider weighted redshift distributions, and use 
the mean of these weighted distributions to give a represen- 
tative redshift for each slice. Explicitly, 



ZA = 



(16) 



(14) 



and similarly for zb- The values for za, zb and the boundary, 
z spin are given in Table [2731 



3 MEASURING THE CLUSTERING 

Power-spectra for the galaxy samples were c alcula ted us- 
ing the method described by iFeldman et al.l l| 19941 ). Each 
galaxy distribution is converted to an over-density field by 
placing the galaxies on a grid and subtracting an unclustered 
"random catalogue" , which matches the galaxy selection. To 
calculate this random catalogue, we model the redshift dis- 
tribu tion of the galaxies using a spline model (|Press et al.l 
1992), and mode l the angular mask using a routine based 
on a HEALPIX dGorski et al.ll2005 | ) eq ual-area pixelization 
of the sphere as in lPercival et al.l (2009). 

Each galaxy and random is weighted using the method 
described in Section 12.2.31 In order to assign luminosities 
to the randoms, we randomly draw a luminosity from the 
galaxy catalogue, with the constraint that the redshift of 
the data point lies within 0.02 of the redshift of the new 
random point. The weight of the same galaxy is assigned to 
the random point - this ensures that the redshift distribu- 
tions of the Vmaz-weighted samples are the same for the data 
and random catalogues. For simplicity, we do not include a 
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JW r,0.1 


ZA 




Zsplit 


r 0.1 


oampt 
Va 


: size 
Vb 


n - m ci f c ri r>n 

11111<_Hj1_.1H_.1_1. 


-23 


0.44 


0.29 


0.38 


-23.04 


12 998 


11 938 


^-matched 


0.44 


0.29 


0.38 


-23.04 


12,998 


11,615 


n - m ci 1" r" n r>n 

ll lll<_tH_.lH_.l_l. 


-23.9 


0.43 


0.28 


0.37 


-22.95 


19 567 


16 226 


£- matched 


0.43 


0.28 


0.37 


-22.96 


19,567 


15,650 


T"l IT) LI *f f"' rl Or 1 
11 lIldH_.lll_.Ll. 


-23.8 


0.43 


0.28 


0.36 


-22.89 


27 328 


19 433 


^-matched 


0.43 


0.28 


0.36 


-22.90 


27,328 


18,300 


n-matched 


-22.7 


0.42 


0.27 


0.35 


-22.81 


35,632 


23,654 


^-matched 


0.42 


0.27 


0.35 


-22.83 


35,632 


22,135 


n-matched 


-22.6 


0.41 


0.26 


0.34 


-22.75 


43,828 


27,269 


^-matched 


0.41 


0.26 


0.34 


-22.77 


43,828 


25,209 


n-matched 


-22.5 


0.40 


0.26 


0.33 


-22.70 


50,994 


29,384 


^-matched 


0.40 


0.26 


0.33 


-22.73 


50,994 


26,677 


n-matched 


-22.4 


0.40 


0.26 


0.33 


-22.68 


56,632 


29,687 


^-matched 


0.40 


0.26 


0.33 


-22.72 


56,632 


26,155 


n-matched 


-22.3 


0.39 


0.25 


0.32 


-22.65 


60,828 


30,855 


^-matched 


0.39 


0.25 


0.32 


-22.70 


60,828 


26,915 



Table 1. Basic parameters for the LRG samples analysed in this paper. 



luminosity-dependent bias model that normalizes the fluc- 
tuations to th e amp litude of L» galaxies as advocated by 
iPercival et all i|2004h . 

Inclusion of the standard IFeldman et all (|1994| ) weight 

(17) 



l + n>P 

would potentially change the match between the high-z and 
low-z samples, because the distribution of galaxy number 
densities in each slice will be correlated with galaxy prop- 
erties, and the correlation may be different for the high-z 
and low-z samples. We therefore do not include this weight. 
Given that the number density of the total SDSS LRG pop- 
ulation is approximately constant with redshift, this only 
has a small impact on the error with which we can measure 
the galaxy clustering strength. 

Power spectra were calculated using a 1024 3 grid m a 
series of cubic boxes. A box of length 4000 h -1 Mpc was used 
initially, but we then sequentially divide the box length in 
half and apply periodic boundary conditions to map galaxies 
that lie outside the box. For each box and power spectrum 
calculation, we include modes that lie between 1/4 and 1/2 
the Nyquist freq uency (similar to the method described by 
ICole et al.ll2005l l. and correct for the smoothing effect of 
the cloud -in-cell assignment used to locate galaxies on the 
grid (e.g. iHocknev k, Eastwood! Il98ll . chap. 5). The power 
spectrum is then spherically averaged, leaving an estimate 
of the "reds hift-space" power. Th i s method is the s ame as 
that used in IPercival et alj l|2009h : iReid et all {200$ ). 

In order to interpret the evolution in the power spec- 
tra measured between high-z and low-z samples, we need 
to know the effect of the survey geometry on the recovered 
power. This window function can be expressed as a matrix 
relating the true power spectrum evaluated at wavenumbers 
k n , Ptrue(kn), to the observed power spectrum P b s ,i, where 
i refers to the different band-powers with central wavenum- 
bers 

Pabs.i = W ^ k n )Ptrue{k n ) - W(ki, 0). (18) 

n 

The term W(ki,0) arises because we estimate the average 



galaxy density from the sample, and is related to the integra l 
constraint in the correlation function l|Percival et al.ll2007T ). 
The window function allows for the mode-coupling induced 
by the survey geometry. Windo w functions for the mea- 
sured power spectrum (Eq n. 15 of[P ercival ct al. 2004 ) were 
calcul ated as described in IPercival et al.l l|200ll ). I Cole et all 
|2005l ). and IPercival et al l (|2007f T 

The covariance matrices for the measured power spec- 
trum band-powers (including correlations between the power 
spectra from the different redshift sli ces), were calculated 
from 10 4 Log-No rmal (LN) catalogues l|Coles fc Joneslll99ll ; 
ICole et all 120051 '). Catalogues were calculate d on a (512) 3 
grid w ith box length 4000 ft -1 Mpc as in IPercival et al.l 
(2009), where LN catalogues were similarly used to estimate 
covariance matrices. Unlike iV-body simulations, these mock 
catalogues do not model the growth of structure, but instead 
return a density field with a log-normal distribution, similar 
to that seen in the real data. The window functions for these 
catalogues were matched to that of the halo catalogue. The 
input power spectrum was a ACDM model matched to the 
large-scale shape of the observed power spectra. 

After removing the effect of the window, we simulta- 
neously fitted the amplitude of the low- and high-redshift 
power-spectra amplitudes using the full covariance matrix. 
We are interested in the large scales only - small scales will 
be affected by intra-halo terms and are sensitive to merg- 
ing happening at these scales. The results presented in Sec- 
tion [5]3] are amplitudes fitted only up to fc max = 0.15. Our 
results are robust to changing this scale cut-off within the 
range fc max w 0.1 — 0.2. Even though we use the full co- 
variance matrix when fitting the amplitudes (which takes 
into account correlations between the two redshift slices), we 
find that the power-spectra are almost independent across 
the two slices. The 1-sigma errors shown in Section 15.31 are 
therefore simply calculated assuming this is strictly true. 

The IFeldman et al.l (|l994h methodology assumes that 
galaxies form a Poisson sampling of a Gaussian random field. 
The resulting shot noise, which arises because of this sam- 
pling, is traditionally subtracted from the measured power 
spectra. Weighting galaxy populations at different redshifts 
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by luminosity means that the large-scale clustering strength 
is expected to be unchanged by loss-less mergers. However, 
the predicted shot-noise would be changed by galaxy merg- 
ers: for example, the expected shot-noise from two equal- 
weight galaxies, or one galaxy with twice the weight, are 
different. Consequently, for our luminosity weighted sam- 
ples, we consider a scenario where we fit both the large-scale 
amplitude and the shot-noise, in addition to subtracting the 
standard Poisson shot-noise. 



4 PREDICTING THE CLUSTERING 
AMPLITUDE 

The observed evolution in the large-scale galaxy power- 
spectrum is the combination of the growth of the fluctua- 
tions in the underlying density field, parametrized by the lin- 
ear growth factor D(z), and the evolution of galaxy bias set 
by the peculiar velocity field. We assume that the expected 
evolution of the amplitude of the redshift-space pow er spec- 
trum between the two redshift slices is given by (|Kaiserl 
1 19871 ) 

P ^ z ^ = Stfrf + M §rP^,z A ). (19) 

b A + 3°AJA + s/ A V A 

Here / is the standard logarithmic derivative of the linear 
growth factor with respect to the logarithm of the scale fac- 
tor. For simplicity, we take / ~ f2m Z \ = 0.25 with 
7(2:) ~ 0.557 — 0.0 22:, which is accurate to .3% for this 
cosmological model |Polarski fc Gannouii|[200l l. 

For a passively evolving pop ulation, we assume that the 
evolution in the bias is given bvlFrvl (|l996l ) 

b z = (6 - 1)/D Z + 1, (20) 

where 60 is the bias at zero redshift. In practice, we compute 
&a from the data, and calculate 

b B = + 1. (21) 

L>B 

It is worth noting that the iFrvl (|l996T ) model assumes 
that the distribution of galaxy velocities matches that of the 
mass. Simulations show that this is reas onable for the haloe s 
in which LRGs are expected to reside (I White et al.ll2007T ). 
although it is possible to create models, albeit ones that are 
somewhat contrived, in which this would not hold, such as 
placing LRGs at stationary points in the density field. 

To calculate the grow th factor we use the fitting formu- 
lae of ICarroll et alj (ll992T ) to calculate the correction to the 
growth factor from an Einstein-de Sitter to a ACDM model. 

Fig. [6]shows the expected evolution in the bias, growth 
factor, and overall power-spectrum amplitude for the red- 
shift ranges we explore in this paper. 



5 RESULTS 

We now present results based on using our high and low 
redshift samples to perform the tests of passive evolution 
described in Section [1] on the number and luminosity densi- 
ties, the luminosity function, and the clustering strength. 



5.1 Number and luminosity densities 

Matching samples based on total numbers of galaxies, obvi- 
ously allows us to test how the number density changes be- 
tween the redshift slices. If some proportion of the galaxies 
seen in the high redshift sample are expected to merge be- 
fore present day, then the merger products will be included 
in the low-redshift sample. In this situation, matching num- 
ber density between samples will bring extra galaxies into 
the low-redshift sample. We should also find that the total 
luminosity changes between the samples if no light is lost 
to the intra-cluster medium through the mergers. Under the 
same assumption, matching samples based on total luminos- 
ity will not bring in these extra galaxies, although we will see 
a reduction in the total number of galaxies. Table [2] shows 
our results for a variety of magnitude ranges; all ratios are 
of za/zb- 

Fig. [7] presents the values in bold in Table f2] and shows 
how the change in number and luminosity density changes 
as a function absolute magnitude. The behaviour seen is per- 
fectly consistent with a scenario where a small proportion 
of the LRGs within the sample at high-redshift merge to 
give brighter LRGs at low-redshift - we see a decrease in the 
number density for luminosity-matched samples and an in- 
crease in the luminosity density for number-density matched 
samples. This is not, however, the only explanation. Merging 
could have also happened between galaxies outwith the the 
sample - i.e., galaxies fainter than the magnitude cut merg- 
ing into brighter galaxies between the two redshifts. As we 
will see later, our clustering analysis will help us distinguish 
between these two scenarios. 

It is worth pointing out that the interpretation of Fig. [7] 
and Table f2] is not straight-forward because the redshift 
slices change as a function of magnitude. As we include 
fainter galaxies in the sample, these are predominantly at 
low redshift, and alter the median value that is used to split 
the sample. The effect measured is probably a combination 
of the inclusion of fainter galaxies and the fact that the split 
occurs at lower redshift - the latter is a small effect, but both 
can act to increase the number of mergers. 

To make the interpretation easier, we can translate the 
ratio of galaxy number density into a percentage number 
of mergers per Gyr -1 , following the assumption that differ- 
ences in number density are caused by loss-less mergers. This 
only makes sense in our luminosity-matched samples, where 
the low-redshift sample does not contain extra galaxies. We 
take the effective time interval At as being that between za 
and zb, and calculate 

r N (M < M r ,o.i) = (l-— ) (22) 
\ UA J At 

where tia and tib are weighted number densities as de- 
fined in (|12[) and (|13|l . We show this rate, as a percentage 
and as a function of magnitude as the solid line in Fig. [8] 
the error bars are Poisson errors. Clearly, interpreting this 
plot as a merger rate only makes sense if one assumes that 
LRG-LRG mergers (from within each sample) are the cause 
to the change in the number density. If instead we have 
fainter galaxies merging together to enter the sample then 
the change in numbers can only be interpreted as a more 
general growth rate. 

We can also do a similar analysis for the fractional lumi- 
nosity growth for samples that have been matched to have 
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Figure 6. The evolution of the bias (left), growth factor (middle) and power-spectrum large-scale amplitude (right) for a population of 
passively evolving galaxies. The dashed line assumes an Einstein de-Sitter model, and the solid line includes a correction for ACDM. 
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i 


Af r ,o.i 


< 


-23.0 


ZA 
zb 
Ratio 


4.56e-06 
4.56e-06 


8886 
9082 
0.978±0.0109 


4.98e-06 
4.76e-06 


9705 
9522 


4.56e-06 
4.43e-06 
1.026±0.0133 


8886 
8885 


4.98e-06 
4.63e-06 


9705 
9308 


Mf.,0.1 
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-22.9 


ZA 
ZB 
Ratio 


6.44e-06 
6.44e-06 


11761 
12069 
0.974±0. 00949 


7.52e-06 
6.97e-06 
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6.44e-06 
6.24e-06 
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11761 
11761 


7.52e-06 
6.72e-06 


13553 
12634 


Af r ,0.1 
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-22.8 


ZA 
ZB 
Ratio 


8.38e-06 
8.38e-06 


14426 
15023 
0.960±0. 00888 


1.05e-05 
9.36e-06 


17644 
16626 


8.38e-06 
7.97e-06 
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14426 
14426 


1.05e-05 
8.80e-06 


17644 
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-22.7 


ZA 
ZB 
Ratio 


1.04e-05 
1.04e-05 


17080 
17770 
0.961±0. 00588 


1.42e-05 
1.24e-05 


22423 
20724 


1.04e-05 
9.89e-06 
1.052±0. 00826 


17080 
17080 


1.42e-05 
1.156e-05 


22423 
19626 


Mr,o.i 


< 


-22.6 


ZA 
ZB 
Ratio 


1.22e-05 
1.22e-05 


19253 
20102 
0.958±0. 00566 


1.83e-05 
1.54e-05 


27287 
24579 


1.22e-05 
1.16e-05 
1.058±0. 00790 


19253 
19253 


1.83e-05 
1.42e-05 


27286 
23033 


M r , .i 


< 


-22.5 


ZA 
ZB 
Ratio 


1.38e-05 
1.378e-05 


20975 
22105 
0.949±0. 00557 


2.23e-05 
1.78e-05 


31534 
27669 


1.38e-05 
1.28c-05 
1.072±0. 00783 


20975 
20975 


2.23e-05 
1.61e-05 


31534 
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< 


-22.4 


ZA 
ZB 
Ratio 


1.48e-05 
1.48e-05 


21983 
23523 
0.934±0. 00557 


2.54e-05 
1.949e-05 


34543 
29673 


1.48c-05 
1.36e-05 
1.094±0. 00801 


21983 
21983 


2.54e-05 
1.70c-05 


34543 
26757 


M r ,o.i 


< 


-22.3 


ZA 
ZB 
Ratio 


1.57e-05 
1.57e-05 


22811 
24478 
0.932±0. 00552 


2.81e-05 
2.10e-05 


37014 
31501 


1.57e-05 
1.42e-05 
1.099±0. 00798 


22811 
22811 


2.81e-05 
1.81c-05 


370134 
28113 



Table 2. The number and luminosity densities of high (za) and low (zs) redshift samples for a range of magnitude cuts, for different 
matching schemes and using different weights. The values shown for a Vmax weight are given only for reference, and all results are based 
on the V ma tch weights. Magnitudes are K+e corrected to z = 0.1; number densities are in units of Mpc — 3 and luminosity densities are in 
arbitrary units. Note that za and zb are not necessarily the same for the different magnitude cuts. A plot of the ratios in bold in these 
table is shown in Fig. [7] 



the same number density. In this case, we compute 



r e (M < M r ,o.i) = 




where I a and £b are weighted luminosity densities as defined 
in f| 10 p and We show this fractional growth rate as the 



dashed line in Fig. [8] This gain in luminosity to low redshift, 
depends on the new galaxies brought into the low-redshift 
sample, and cannot therefore be as easily interpreted by a 
merger model as the luminosity-matched samples. We quote 
the values of rjv and re, including Poisson error bars, in 
Table El 



LRG evolution in DR7 11 




nqo r ... i ... i ... i ... i . 

-23.2 -23.0 -22.8 -22.6 -22.4 



Figure 7. The ratio of the number density at high and low red- 
shift when matching using the luminosity density (black solid line) 
and the ratio of the luminosity density at high and low redshift 
when matching using the number density (black dashed line). 
This plot shows the values in bold in Table [2] The blue and red 
lines show the estimated effect of uncertainties in the IMF slope 
(see Section \E\ for details) 




Figure 8. The solid black line shows rjv x 100: the percentile 
growth rate in number of galaxies per Gyr, as defined by equation 
1221 This can be interpreted as a merger rate if we assume that 
galaxies in the sample at high redshift have merged. The dashed 
black line shows re X 100: the percentile increase in luminosity 
density for samples matched to have the same number density, as 
defined by equation 1231 Both quantities are shows as a function 
of the faint magnitude limit, and given in Table [3] The blue and 
red lines show the estimated effect of uncertainties in the IMF 
slope (see Section [6] for details). 



5.2 Luminosity function 

The next natural step is to construct a luminosity function 
and study its evolution, which we compute for each of red- 
shift slices as 

0(M r ,o.i)AM P| o.i = TT~- (24) 

Mi6AM r , .i shce 

This does not give a luminosity function in the traditional 
sense, but rather it gives a population-weighted luminosity 



1V1 r,0.1 


rjv x 100 


re x 100 


-23.0 


2.04 ±1.035 


1.72 ±0.847 


-22.9 


2.43 ±0.883 


2.02 ±0.735 


-22.8 


3.82 ±0.812 


3.20 ±0.689 


-22.7 


3.81 ±0.640 


3.13 ±0.455 


-22.6 


4.24 ±0.614 


3.43 ±0.441 


-22.5 


5.23 ±0.612 


4.21 ±0.436 


-22.4 


6.74 ±0.630 


5.51±0.438 


-22.3 


7.10 ±0.628 


5.75 ±0.435 



Table 3. The measured values of rjv and re, as defined by equa- 
tions J22]) and J23]l and plotted in Fig|8] 
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Figure 9. The luminosity function in Va (solid line) and Vb 
(dashed line), as defined by Equation I24I The redshift boundary 
between the two volumes is z = 0.33. Given our V^at^ weight- 
ing scheme, this luminosity function is not corrected for volume 
completeness. Instead, it is constructed such that the population 
at low redshift is perfectly matched to the population at high 
redshift, see text for details. 



function. I.e., the luminosity function of galaxy samples at 
two redshifts that represent the same population of galaxies 
in equal terms, albeit perhaps incompletely. We show this 
pseudo-luminosity function in Fig. [9] 

Fig. [9] shows that the information content is not sig- 
nificantly improved from that in Figs. [7] & [5] which show 
the change in the full sample. This is simple to understand: 
small differences in number density are always diluted when 
one splits them into magnitude bins. Fig. [9] does make it 
clear that one needs a way to enrich the low-redshift sample 
with galaxies of M r> o.i > —23, by means other than passive 
stellar evolution. 



5.3 Clustering 

The last test that we perform is based on measuring the evo- 
lution of the clustering of our matched samples of LRGs. The 
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Vinatch -weighting makes sure that we are comparing like- 
with-like galaxies for clustering measurements. Matching 
the total luminosity means that, even allowing for loss-less 
mergers, we are comparing galaxies at low-redshift that are 
evolved products of the high-redshift galaxies without bring- 
ing extra galaxies into the sample. Additionally, weighting 
by evolved luminosity, means that the large-scale clustering 
should not be affected by these mergers: on large-scales there 
is no change in the clustering strength if two galaxies merge. 
The method for power spectrum measurement was described 
in Section [3] We compare this carefully constructed test 
against the clustering evolution observed for the samples 
matched by number density. The passively evolving model 
that we test against was described in Section [4] 

We fit an amplitude to the large-scale power-spectra at 
each redshift, as described in Section[3] We are interested in 
the evolution of this amplitude from high- to low-redshift, 
so we will plot ratios of these amplitudes. Fig. [IJJshows this 
ratio as a function of limiting magnitude at high-redshift, for 
the luminosity-matched and number density matched sam- 
ples as the black solid line and stars. 

There are four distinct points to appreciate about 
Fig. [10] 

(i) neither the luminosity-weighted nor the standard- 
weighted power-spectrum follow passive evolution at the 
faint end; 

(ii) in both cases the most luminous objects are found to 
be consistent with passive evolution; 

(iii) the departure from passive evolution is comparable 
in both cases; and 

(iv) in both cases, there is an under-evolution of the 
power-spectrum amplitude with redshift, with respect to 
that expected from passive evolution. 

The departure from passive evolution is less clear in the 
luminosity-weighted power-spectrum when we fit for shot- 
noise independently at high- and low-redshift (in the open 
black squares), which we do for the reasons given in Sec- 
tion 21 The general behaviour is, however, maintained. A 
clearer signal is at the moment not possible due to the lim- 
ited redshift baseline given by our sample. 



6 MODEL DEPENDENCE 

Our results are naturally dependent on the choice of stellar 
model that describes the passive evolution of stellar colours 
and luminosity. We justify our choice of model by noting that 
it describes the colour evolution of at least a sub-sample of 
LRGs (those tested in iMaraston et alJboOSl ) over a redshift 
range that goes beyond what is needed here. Implicitly, we 
assume that this model is a good description al all LRGs 
in our sample. In this section we explicitly discuss possible 
ways in which the modelling may be wrong. 

The luminosity evolution, for one, is not as well con- 
strained as the colour evolution. We must keep in mind that 
a model that would predict a different rate of stellar fad- 
ing would impact on our weights that in turn would impact 
on the sample of galaxies that is returned by our matching 
scheme. The dominant factor in this case is t he slope of the 
Initia l Mass Function (IMF) at around IMq (|Conrov et al.l 
2009), but there are in principle other reasons for a change 



in the luminosity evolution. Although quantitatively moti- 
vated by changes to the IMF, the test we present in this 

section is general for all. 

As explicitly demonstrated in Fig. 8 of IConrov et al.l 

(2009), a change in the slope of the IMF of around ±0.8 
affects the expected B magnitude change of passively evolv- 
ing galaxy by as much as 0.4 magnitudes for z < 1. The 
details of this calculation most certainly depend on the ex- 
act modelling of the passively evolving galaxy, as well as on 
the studied spectral region. However, we can certainly test 
for an error of this order of magnitude in our samples. We 
assume that the quoted uncertainty is uniformly distributed 
since z = 1 and that \dM\ = 0.2dz. We therefore change 
Equation [TJ to 

M r , . 1 = r p -5log 10 ^^^-^ -Ke(z,0.1)±0.2z (25) 

and re-run our analysis as before. For our redshift slices, 
which have approximately dz — 0.15, the introduced mag- 
nitude difference between high and low-redshift will be of 
the order of ±0.03 magnitudes with respect to the assumed 
IMF slope. Note that what is important here is the differen- 
tial effect with redshift, which affects galaxies in Va and Vb 
with different degrees - we are insensitive to an overall shift 
in luminosity. In the case of a negative shift (dM = — 0.2dz), 
galaxies at lower redshifts will have faded more with respect 
to galaxies at high-redshift, and we will need more objects 
in Vb to explain the luminosity in Va- So whereas before we 
saw a deficit in the number density when matching the light, 
we should now find more closely matched number densities. 
The same reasoning applied to a positive shift (dM = 0.2z), 
means we should require fewer low redshift galaxies. Note 
that when we match by number density the main difference 
comes not from selecting different galaxies, but rather from 
a different evolution of the magnitudes of these galaxies. The 
reasoning is similar to the above - for a negative shift we ex- 
pect the luminosity density at low redshift to be decreased 
with the respect to the fiducial model. So when we saw an 
excess in luminosity when matching by number density, we 
now expect this excess to be reduced. 

We show the effect of a ±0.2z change in the fiducial lu- 
minosity evolution on the ratios of number and luminosity 
densities, as well as growth rates, in Figs [7] and [5] In each 
case, the red (dM = 0.2z) and blue (dM = — 0.2jz) lines 
show the expected evolution of number and luminosity den- 
sities if the luminosity evolution was off by ±0.2z. One can 
see that if our description for the luminosity evolution was 
off by —0.2z, then we could just about reconcile the obser- 
vations with the passive evolution scenario (except perhaps 
for our faintest sample). In other words, we would be en- 
riching the low-redshift sample with more (lower-luminosity) 
objects. However, the galaxies removed or included by ap- 
plying this evolutionary change have the lowest luminosities 
and are the least biased. This change will consequently also 
affect the amplitude of the clustering signal.. We show the 
power-spectrum evolution results in Fig. If 01 The enrichment 
of Vb with lower-luminosity objects causes the amplitude of 
the clustering to be further reduced, as one would expect. 
Therefore, even though a potential error in the modelling 
of the luminosity evolution could perhaps explain the evolu- 
tion of number and luminosity densities without departing 
from dynamical passive evolution, it could not explain the 
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Figure 10. The ratio of the large-scale amplitude of the power-spectra at high- and low-redshift, as a function of the limiting magnitude 
at high-redshift. Left panel: luminosity-weighted samples and power-spectrum. Right panel: number density samples and uniformly- 
weighted power-spectrum. In both panels the dotted line shows the expected evolution from a pure passively evolving population of 
galaxies. The open squares in the left panel (slightly offset in magnitude for clarity) show the effect of fitting the shot noise separately, 
using a k max = 0.2 - see Section [3] for details. The blue and red lines show the estimated effect of uncertainties in the IMF slope (see 
Section [6] for details). These data points are highly correlated, as each fainter sample includes the brighter ones. 



evolution of the clustering. It is worth commenting on two 
further aspects of Figs [7] [8] and 1101 In all cases, we see a 
stronger effect in the case of a — 0.2z shift, and one which 
increases to fainter magnitudes. This is probably explained 
by the slope in the luminosity function at low-redshift. By 
adding —0.2z, we are adding objects to the low-z sample 
(whilst keeping the high-z sample more or less the same - 
see next paragraph), and the results will depend heavily on 
the luminosity and number densities of these objects (which 
evolves faster towards fainter magnitudes). The other aspect 
worth mentioning is that we do not necessarily expect the 
blue and red lines to bound the solid black line in Fig. 1101 
In this case, the change in the results is given exclusively by 
the properties of the objects that either enter or leave the 
sample (with respect to the black line), and given that their 
number is small then this is effectively a noisy measurement. 
The exception is given by the increasing offset of the blue 
line, for which this number is in fact larger for the reasons 
mentioned above. 

There is a final subtlety in the above analysis. When we 
apply Equation (|25p . we are effectively changing the abso- 
lute magnitude of all the objects in our sample. If we apply 
the same magnitude cuts as before, we then select a differ- 
ent sample of galaxies and comparison proves difficult. We 
therefore choose to shift the absolute magnitude limits at 
high-redshift by the same amount, in order to keep roughly 
the same objects in Va for all cases. Note that both ap- 
proaches would be valid - in one case we would be studying 
the effect of Equation (|25|) for fixed magnitude bins, and in 
the other we are studying the effect of Equation (125[l on a 
fixed sample of galaxies. 

We also do not directly address the question of contam- 
ination. As we find a clear deviation from passive evolution 
in at least part of our sample, we should keep in mind that 
this could very well be because the galaxies that reside in the 
lower-mass haloes are indeed not LRGs in the traditional, 
stellar population, sense of the word - i.e., their colours may 



be mimicking those of LRGs due to dust, for example. This 
is perfectly OK - we are test ing whether the galaxi es that 
are selected according to the lEisenstein et al.l (|200ll ) target 
selection algorithm are passively evolving in a dynamical 
sense, independently of what causes their colour evolution. 
We must only worry that we have a good model for the 
description of its colour evolution with redshift. 

Finally, let us also consider a scenario in which LRGs 
are generally passively evolving, but their formation epoch 
has a dependence on luminosity, or mass. This scenario is 
easily motivated by the literature, as several authors have 
found a dependen ce of mean age with luminosity in early- 
type galaxies (e.g. ICaldwell et al.ll2003l ; iThomas et ai1l2005l ; 
ICiemens et al.l 12006). and would effectively mean that the 
M09 model, even if correct for the most massive LRGs, may 
become increasingly inadequate for fainter objects. However 
even the LRGs in our faintest sample, with M r< o.i > —22.3, 
have typical masses of 1O 11,5 M0, and velocity dispersions 
greater than 200km s . Given these numbers, we should 
expect a difference in formation age of less than lGyr, which 
would make little difference for the spectral evolution at the 
redshifts we are considering. This gives us confidence that 
the interpretation we present in the next section is driven by 
the dynamical evolution of the objects in our sample, rather 
than by an increasingly inadequate model as we go to fainter 
galaxies. 

The dependence on the stellar model is characteristic of 
all studies that need to match LRG samples at low and high- 
redshift, and emphasises the importance of having a good 
model for the colour and luminosity evolution of LRGs. 



7 INTERPRETATION 

In Sections l5.1U5.2l and l5.3l we have shown the results of test- 
ing passive evolution of LRGs using number and luminosity 
density evolution, as well as the evolution of the clustering 
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of LRGs. If our goal is to simply test dynamical passive evo- 
lution of LRGs, then our three measurements give a very 
clear answer - LRGs, as selected, do not all follow passive 
evolution. Fig. [71 [8l [9l and [lOl consistently show evidence for 
merging of some sort. In all cases we also see a clear de- 
pendence on luminosity - the brightest galaxies show the 
smallest departure from pure passive evolution. 

To explain the change in weighted number and lumi- 
nosity densities seen in Fig. [7J and Fig. [5] we have to call in 
on some kind of merging scenario that explains the increas- 
ing luminosity density for samples with the same number 
density, or the decrease in the number of galaxies needed to 
explain the same amount of light. This can happen in three 
distinct ways: 

(i) merging within the sample; 

(ii) galaxies outside the sample merging with galaxies in 
the sample; or 

(iii) galaxies outside the sample merging together and 
getting bright enough to make it into the sample. 

The clustering result is the least clear, because of the 
relatively small redshift range probed. At the same time, 
however, the introduction of the luminosity-weighted power- 
spectrum has the most potential to tell us something about 
the nature of the evolution of LRGs, and differentiate be- 
tween the three scenarios above. 

If we restrict ourselves to the samples matched by 
number density and a uniformly-weighted power-spectrum, 
then our results ar e in agreement with previous work from 
IWhite et all (120071 ); IWake et al.l (120081 ). where it was found 
that the evolution of the large-scale bias or clustering am- 
plitude show very little evolution with redshift. A valid in- 
terpretation for this signal is the merging of objects at high- 
redshift into one LRG at low-redshift, which decreases the 
number density of objects in high-mass haloes and brings 
the large-scale clustering amplitude down. 

However, if the evolution (or lack of) seen in the right- 
hand panel of Fig. [10] is explained by the merging of satellite 
LRGs within the sample, then we expect the luminosity- 
weighted power-spectrum to follow passive evolution. The 
fact that the data in the left-hand panel matches that in 
the right, showing some evidence for a departure from this 
hypothesis, leads us to consider a scenario where the lack 
of evolution in the large-scale power is due to new objects 
entering the sample between the two redshifts. 

We can further assert how the new objects are enter- 
ing the low-redshift sample. If the growth happens at the 
massive end, we expect the evolution of the power-spectrum 
to overshoot that expected from passive evolution, as more 
weight is given to objects that are intrinsically more strongly 
clustered. In reality, we observe the opposite behaviour. The 
only explanation that is consistent with both measurements 
is that objects in lower mass haloes are gaining weight from 
high to low redshift. 

We cannot distinguish between a scenario where galax- 
ies that fall completely outside the sample at high-redshift 
come together and become bright enough to make it into 
the sample at low redshift, from a scenario where the faint 
end of the sample at high redshift is gaining weight from the 
merging of companions that initially fall outwith the sam- 
ple. These two scenarios, however, are one and the same - 
the growth of LRGs, as a population, is happening at the low 



mass end. As we slide our magnitude limits we probe differ- 
ent regions of this growth but we observe it increases as we 
go down in halo mass/luminosity. 



7.1 Intra-Cluster Light 

We have so far ignored the possibility that light is lost in a 
merging event. If this happens then the luminosity- weighted 
power-spectrum is no longer expected to follow passive evo- 
lution in the case of merging within the sample. We can ask 
the question: is there a fraction of light loss per merger of 
given mass that can explain the observed lack of evolution 
in the luminosity-weighted power-spectrum? To answer this 
question quantitatively, we need to know how the amplitude 
of the power-spectrum changes as a function of luminosity, 
P(k, L), which we leave for a follow-up paper. 

Qualitatively however, we can still make the following 
observation: we need to increase the weight of objects resid- 
ing in smaller haloes at low redshift to explain the observed 
deficit in power. It follows that any light loss at the bright 
end would have to be such that the resulting increase in 
luminosity of the bright objects is low enough that the clus- 
tering of these objects does not dominate the overall cluster- 
ing signal. This would point towards a differential mass loss 
fraction with luminosity, with the most luminous objects 
losing the most mass to the ICM. The slope of this relation 
is related to the slope of dP(k, L) /dL, but we would still 
require preferential mass growth at the low-mass end - in 
fact, even more so. It therefore seems inevitable to conclude 
that the LRG growth is happening predominantly in lower 
mass haloes. 



8 COMPARISON WITH PREVIOUS WORK 

Our analysis has given us as a behavioural description of 
the growth of LRGs. Under our hypothesis, the luminos- 
ity growth that we present in Fig. [5] is being introduced 
by objects that live in smaller mass haloes, rather than 
satellite accretion into luminous central objects. This means 
that we cannot interpret the changes in number density be- 
tween high and low redshift as a merger rate. The luminosity 
growth, in turn, is more directly comparable with previous 
results. 

Even so, the comparison with other work is far from 
straightforward, given the different colour/luminosity selec- 
tion, number density, redshift range and absolute magni- 
tudes explored in each work. Most simply, we expect LRG 
growth to increase with redshift, and with decreasing lumi- 
nosity. Table [5] presents a non-exhaustive collection of recent 
literature results that measured either the merger rate, or 
the luminosity growth rate of LRGs using a variety of tech- 
niques. The results that we present in Fig. [H] can at least be 
said to sample the same range of valu es in Table El with th e 
exception of the growth measured bv lMasiedi etalJ (2008). 
This may be reconciled with our interpretation if we con- 
sider that we need new objects becoming bright enough to 
become LRGs form high to low redshift and that is what 
dominates LRG growth. By cross-correlation LRGs with the 
main galaxy population, the authors concentrate only on the 
growth of existing LRGs. 
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Publication 



Redshift range 



Luminosity Growth 
(Gyr- 1 ) 



Merger rate 
(Gyr" 1 ) 



Number density 
(Mpc" 3 ) 



Masiedi et al. (2006) 0.16 - 0.36 

Wake et al. (20061 0.2 - 0.55 

Brown et al. (2007) < 0.9 

White et al. (20071 0.5-0.7 

Cool et al. (20081 0.1 - 0.9 

Masiedi et al. (20081 0.16 - 0.30 

Wake et al. (20081 0.2 - 0.55 
De Propris et al. (20101 0.45 - 0.65 



«3% 

< 18% 
6.8% 

< 1.7/i% 



0.6 x 10 4 Gpc" 3 
< 10 % 

a 3.4% 



2.4 % 

0.8 x 10 4 Gpc" 3 



< 1-25 x 10~ 6 
; 3 x 10~ 4 
; 10~ 3 /i 3 
; 3 - 9 x 10~ 5 

; 3 x 10~ 4 



Table 4. A summary of the values obtained in the literature for the growth and merger rates of LRGs. These may be compared with 
the results we present in Tableland in Fig|8] but note the varying redshifts and number densities in each study. 



The luminosity-weighted power-spectrum is a power- 
ful tool to disentangle merging scenarios: mergers between 
satellites and centrals without light loss would give differ- 
ent evolution in the large-scale clustering strength between 
number density and luminosity matched catalogues. Instead 
we see similar evolution. Additionally, in order to match the 
observed decrement in the power spectrum evolution for the 
less-luminous galaxies and the increased rate of evolution in 
number and luminosity density, the simplest explanation is 
that LRGs need to be introduced in smaller haloes which 
are intrinsically less clustered. This interpretation is at odds 
with the one typically offered by the halo model, which is 
based on the small-scale clustering of samples not split by 
luminosity and that are matched by number density. In the 
"standard" halo-model based explanation, departures from 
passive evolution are due to satellite-central mergers. The 
need to match by number density when applying the halo 
model to test evolution is on itself a problem in performing 
these fits as, if the model suggests some fraction of satellites 
merge onto the central galaxy, then the number density at 
low redshift must also be reduced. How exactly to do this is 
not immediately obvious, and has been done differently by 
different authors. 



9 DISCUSSION AND SUMMARY 

We have, for the first time, applied luminosity weighting to 
sample selection and to the power-spectrum, and have pre- 
sented a new method with which to interpret the evolution 
of LRGs. We have also introduced a new weighting scheme 
that allows us to keep most of the galaxies in the sample. 
This has allowed us to study the evolution of LRGs as a 
function of luminosity with unprecedented resolution. We 
have done this by measuring evolution of the number and 
luminosity density of SDSS LRGs with 0.15 < z < 0.5, as 
well as the evolution of their clustering. 

We can summarise our interpretation and conclusions 
in the following bullet points: 

• The evolution of LRGs, as a population, is inconsistent 
with passive evolution. 

• Departure from pure passive evolution is strongest for 
fainter LRGs; bright LRGs are consistent with pure passive 
evolution. 

• We see a lack of evolution in the large-scale luminosity- 
weighted power spectrum for objects with M r ,o.i < —22.8, 



relative to what is expected from passive evolution. This 
effectively rules out option (i) in Section [7J as an explana- 
tion for the merger and luminosity growth rates presented 
in Table [31 although the evidence for this interpretation is 
not strong as a consequence of the relatively narrow redshift 
range probed. 

• To explain this lack of evolution instead we propose that 
LRGs in smaller mass haloes must gain weight (luminosity) 
since z — 0.5. It is unclear whether options (ii) or (iii) of 
Section [7] are dominant, but given our sliding magnitude 
cuts the two are effectively the same process. 

• Our interpretation relies on the assumption that light 
is conserved when two galaxies merge. However, even if this 
is not true, any weight given to LRGs in high-mass haloes 
must be offset by new objects in low-mass haloes entering 
the sample at z < 0.5. This, however, only increases the 
need to introduce less clustered objects into the low-redshift 
sample. 

• Our results are not inconsistent with the halo model 
interpretation, nor previous work, if we restrict ourselves to 
the same observables. The added information in this paper 
comes from the matching of the samples on luminosity den- 
sity (rather than number density) and on the evolution of 
the large-scale luminosity-weighted power-spectrum. 

• We have explicitly estimated the effect of uncertain- 
ties in the IMF slope in the rate of fading of the stellar 
populations and in the subsequent selection and interpreta- 
tion of our samples. The departure from passive evolution 
seen in the number and luminosity densities can certainly 
be attributed to an (extreme) uncertainty of the IMF slope. 
However, the resulting clustering signal is then clearly in- 
consistent with dynamical passive evolution. 



In terms of pushing this analysis fu rther, the Baryon Os - 
cillation Spectroscopic Survey (BOSS: ISchlegel et al.| [2009). 
part of the SDSS-III project, will measure redshifts of LRGs 
out to z < 0.7. This will provide the lever-arm required to 
fully test the evolution of clustering strength as a function 
of redshift. When combined with luminosity- weighting this 
will allow us to test whether our interpretation is indeed 
correct. 
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